Quantum resources for hybrid communication via qubit-oscillator states 



Tommaso Tufarelli 1 , Davide Girolami 2 , Ruggero Vasile 2 , Sougato Bose 3 and Gerardo Adesso 2 
1 QOLS, Blackett Laboratory, Imperial College London, London SW7 2BW, United Kingdom 
2 School of Mathematical Sciences, The University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom 
3 Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom 

We investigate a family of qubit-oscillator states as resources for hybrid quantum communication. They result 
from a mechanism of qubit-controlled displacement on the oscillator. For large displacements, we obtain ana- 
lytical formulas for entanglement and other nonclassical correlations, such as entropic and geometric discord, 
in those states. We design two protocols for quantum communication using the considered resource states, a 
hybrid teleportation and a hybrid remote state preparation. The latter, in its standard formulation, is shown to 
have a performance limited by the initial mixedness of the oscillator, echoing the behaviour of the geometric 
discord. If one includes a further optimization over non-unitary correcting operations performed by the receiver, 
the performance is improved to match that of teleportation, which is directly linked to the amount of entangle- 
ment. Both protocols can then approach perfect efficiency even if the oscillator is originally highly thermal. We 
discuss the critical implications of these findings for the interpretation of general quantum correlations. 
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I. INTRODUCTION 

Understanding, identifying and exploiting "quantumness" 
in composite systems represent essential steps to grasp the 
fundamental implications of quantum theory Q, and are of 
particular relevance in the race for efficient information and 
communication technology applications defying the classical 
boundaries [2]. Nonlocality and entanglement are clearcut 
signatures of nonclassicality OIU. However, it has recently 
been acknowledged that, focusing on aspects of correlations 
among quantum systems, manifestations of quantumness can 
exist even in absence of entanglement [5]. Quantum discord 
and related measures of general quantum correlations are re- 
ceiving widespread attention [ 6 ] as they promise to enable 
a supraclassical speedup in computational frameworks where 
entanglement is not robust enough to endure |7|. Some pro- 
tocols have been identified that appear to take advantage of 
discord-like correlations, rather than entanglement, for their 
functionality (8) ; still, basic questions about the interpretation 
of discord remain unanswered. 

Some steps have been undertaken to understand the inter- 
play between entanglement and general quantum correlations 
ID |9). A hierarchic relation is proven for two qubits |[TQll , 
involving the negativity N ifTTTl as an entanglement measure, 
and the so-called geometric discord T>q lfl2l as a nonclassical- 
ity indicator: V^g > N. Less is known about the structure 
of general quantum correlations in high-dimensional systems, 
with the exception of continuous variable Gaussian states lTT3l — 
16]. 

Here we study hybrid bipartite systems comprising a two- 
level system A (qubit) and a harmonic oscillator B (qumode). 
We consider a class of states pab in which correlations are 
induced by the action of qubit-controlled displacements on the 
oscillator [17]. These states can be implemented in several 
setups fT8l and the required interaction can be exploited for 
universal quantum computation lfT9l . 

The main aim of this paper is construct protocols for quan- 
tum communication using these resource states, and analyze 
their performance in connection with the contents of different 




FIG. 1: (Color online). Schemes for hybrid teleportation (a), and 
hybrid remote state preparation (b), via the shared resource states 
p AB of a qubit A and an oscillator B. Details are provided in the text. 



types of correlations in the states pab- 

We define two state-transfer protocols based on the shared 
states pab, see Fig. [T] One is a hybrid teleportation scheme 
l20l where Bob can teleport an unknown qubit state |0fj> to 
Alice, and whose fidelity is proven to approach unity for a 
particular subclass of resource states. The second is a hybrid 
remote state preparation protocol l2Tll where Alice can mea- 
sure the qubit to remotely prepare Bob's oscillator in some 
(known to Alice) state |<A|° al ). In this case, if the receiver Bob 
can only perform unitary corrections, the fidelity is bounded 
in general by the initial purity of the oscillator [22-24]; how- 
ever, if we allow Bob to perform non-unitary corrections, the 
figure of merit can increase and match the one associated to 
teleportation. 

We thus analyze in detail the nature of correlations in the 
states pab and their role for the performance of the differ- 
ent protocols. We obtain analytical formulas for negativity 
N [11], entropic quantum discord T>z [5], and geometric dis- 
cord T>g lfl2l of the states pab- We find that such states can be 
maximally entangled and maximally discordant in the limit of 
large displacements, while their geometric discord is limited 
by the initial purity of the oscillator and can be thus arbitrar- 
ily small. We argue that this is a direct consequence of the 
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particular geometry of the state- space induced by the Hilbert- 
Schmidt norm, which enters the definition of T>q 1 12 ], lead- 
ing us to conclude that T>q cannot be regarded, in general, as 
a proper 'measure' of nonclassical correlations, in agreement 
with the conclusions of Ref. l25l . 

This work provides three main advances. First, from a 
practical viewpoint, it presents workable protocols for hy- 
brid quantum communication, which can be useful as build- 
ing blocks in any light-matter interfaced implementation of 
quantum information processing [19, 26, 27 ]. Second, from a 
technical viewpoint, it provides useful methods for the analyt- 
ical evaluation of correlation quantifiers in bipartite systems 
with subsystems of different dimensionality l28l . Third, from 
a physical perspective, it exposes the need for a mathemati- 
cally sound and physically meaningful approach to delve into 
the nature and structure of general nonclassical correlations in 
quantum states. 

The organization of the paper is as follows. In Sec.|n|we de- 
fine the studied model and calculate analytically various mea- 
sures of nonclassical correlations for the states pab- In Sec. [111 
we define protocols for teleportation and remote state prepa- 
ration using the shared resource states pab and calculate their 
fidelity, relating it to the correlations present in the states. In 
Sec.|lV]we draw our conclusions. Some technical derivations 
are deferred to Appendixes. 



II. THE STATES AND THEIR CORRELATIONS. 

We consider a family of qubit-oscillator states pab = 
Uab(J3)(Pa ® P°b) U ab& obtained by applying the unitary 
Uab(J3) = Dip-^P) to an initially uncorrected general state 
Pa®P°b- ^ ms m t erac ti° n induces a qubit-controlled displace- 
ment on the oscillator, that is, a displacement where the sign 
of the parameter +J3 is determined by the eigenvalue of the 
qubit Pauli operator cr 3 (17] [T8l . This corresponds e.g. to 
the evolution of the joint system via a coupling Hamiltonian 
H oc cr^ib + M), which can be realized experimentally in a 
number of setups (HI [29). This type of interaction has rele- 
vant applications for the state reconstruction of oscillator net- 
works probed by a single qubit l30l and for quantum compu- 
tation based on light-matter interfaces H91 . We describe the 
state of qubit A before the interaction as = ( ? x r _^j in the 

standard basis {\e), with < p < 1 and |r| 2 < p(l - p). 
The states of the hybrid system after the interaction take the 
form 

Pab = p\e)(e\®D(J3)p B rt(J3) + a-p)\gXg\®rt(0)p° B D(fi) 
+ r\eXg\^D{fi)plD{p) + r\Xe\^D^)plD^). (1) 

We focus on the regime of large displacements — > oo, 
which in practice means large enough such that the overlap 
between the two phase-space domains, associated to p° B dis- 
placed by j3, and to p° B displaced by -/3, becomes negligible. 

Let us calculate the correlations in the states of Eq. ([I]). 
Nonlocality properties have been studied in |3H . 

Entanglement can be quantified by the (normalized) nega- 
tivity CD N(pab) = (\\p\ A B \h - 1), where ||M||i = Tr|M| de- 
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FIG. 2: (Color online). Plot of negativity [Eq.|2) (solid line), lower 
bound on entropic discord [Eq. ([4])] (dashed line), and geometric dis- 
cord [Eq.[5] for jj° b = 0.5, 0.1 (dotted and dot-dashed line, respec- 
tively) calculated for the states p AB with p = 1 /2, as a function of |r|. 
The plotted quantities are dimensionless and correspond to \fi\ — > oo. 



notes the trace norm and p\ B is obtained by partial transposi- 
tion with respect to the qubit only l32ll . We find the following 
rigorous result, whose proof is provided in Appendix |A| 



lim N{pab) = 2\r\ . 



(2) 



Independently of the initial state of the oscillator (namely, 
no matter how thermal it is), one can always find a \/3\ large 
enough such that the qubit-oscillator states are asymptotically 
maximally entangled, when the qubit is initially in a pure 
equatorial state, p = |r| = 1/2. Then the state pab reproduces 
a proper Schrodinger cat state, if one interprets the qubit as 
the 'microscopic' degree of freedom and the qumode as the 
'macroscopic' one l33l . 

Let us now consider more general types of nonclassical cor- 
relations. The geometric discord Dg(Pab) quantifies how far 
(in Hilbert-Schmidt norm) a bipartite quantum state is from 
the set of classical-quantum states lfT2l . It can be measured ex- 
perimentally with direct, non-tomographic methods 134U361 . 
Given a quantum state pab of a (D 2 0(D d system, with A being a 
qubit and B being an arbitrary (finite or infinite) d-dimensional 
system, the normalized geometric discord Dg(Pab) is defined 



as 1121123 D G (p AB ) = 2infn A \\p AB - U a (pab)\\ 2 2 , where the 
infimum is over all von Neumann measurements Ha = {11^} 
on the qubit A with U a (pab) = 1^(1^ Wpab^a ® 1 B ), 
and ||M|| 2 = ^Tr(MM^) is the Hilbert-Schmidt norm. In Ap- 
pendix [B] we provide a useful and compact analytical frame- 
work for the calculation of T>q in 2 x d systems including the 
case d - oo, which we need here. We obtain 
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lim D g (pab) = V B W 



(3) 



The geometric discord can be made arbitrarily small by de- 
creasing the initial purity p° B = tr fi [(p^) 2 ] of the oscillator (see 
Fig. [2]). This result shows that the conjectured ordering rela- 
tion V25g ^ N, which holds for all two-qubit states fTQi is 
violated when at least one subsystem of a bipartite system has 
large dimension, as in our case. For the states pab, the order- 
ing is actually reversed, revealing the quirky situation of states 
with possibly maximum entanglement yet asymptotically van- 
ishing (for p° B — > 0) geometric discord. 

We then ask: are those maximally entangled states infinites- 
imally close to the classical-quantum border in other met- 
rics? To address the question, we evaluate the conventional 
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entropic quantum discord T>z on the states pab- We re- 
call the definition, D z (pab) = M Ua {?(Pab) - I[U A (p AB )]}, 
where I denotes the quantum mutual information. We ob- 
serve that, by means of a local operation on Bob's side, one 
can map pab onto an effective two-qubit state Pab- Entangle- 
ment measures [3 ] as well as the discord T>z with measure- 
ments on A l38l (but crucially not Dg (39)) are monotonic 
under such operations on B, so that by evaluating those corre- 
lation measures on p A B one obtains lower bounds to the cor- 
responding measures for pab ES EH EH- In particular, as 
proven in Appendix [C| Bob can choose two orthonormal vec- 
tors \e), ||) and design a local operation which, in the limit of 
large displacements, converts the state pab to the 'digitalized' 
PaI = P\eeXee\ + r\ee\gg\ + r*\gg\ee\ + (1 - p)\ggXgg\, with- 
out resorting to post-selection. This effective two-qubit state 
can be achieved independently of the initial purity of B. The 
entropic discord in the digitalized state reads 

£>z(PaD = -plog 2 (p)-(l-p)log 2 (l-p) (4) 
+ i[\og 2 (p(l-p)-\r\ 2 )^aog 2 {^)l 

with ( = 7(1 -2p) 2 +4|r| 2 . In the case p = \r\ = 1/2, the en- 
tropic discord D z (p^|) converges to 1 like the entanglement, 
implying Dz(pab) — > 1 in the original states of Eq. ([T]) as well. 

This demonstrates that the Hilbert- Schmidt norm is gener- 
ally unsuitable for defining quantitative distance-based mea- 
sures of correlations — as recognized for entanglement in l42l 
and very recently for geometric discord in (25l — with its de- 
ficiency being even more critically exposed in systems with 
large dimension. We can then reassess some findings in 
the recent literature, such as Dg failing to capture the re- 
source power of the discrete quantum computation with one 
bit [Q21 [36), as evidences of this deficiency. Nevertheless, 
when the purity of the states is fixed, and/or when the dimen- 
sion of the Hilbert space is small enough (e.g. for two qubits), 
the geometric discord returns a reliable quantification of non- 
classical correlations (T0|[24l|43). In general, given its com- 
putability and experimental accessibility (34] [52), it can still 
play a useful role if one correctly regards it not as a measure 
by itself, but as a valid lower bound to regular measures of 
nonclassical correlations like the relative entropy of discord 
(44l[45) (alias one-way deficit (46)). Such a bound becomes 
looser with increasing dimension, as evidenced by our analy- 
sis, and as emerged in the study of continuous variable Gaus- 
sian states (16). 



III. HYBRID QUANTUM COMMUNICATION 
PROTOCOLS 

We now question whether the different aspects of correla- 
tions identified in the states of Eq. ([I]) (see Fig.[2]for a compar- 
ison) can be endowed with operational meanings. We design 
two hybrid quantum communication protocols which employ 
the family of states pab, with Alice operating the qubit A and 
Bob operating the oscillator B, as shared resources. 



A. Hybrid teleportation 

Alice and Bob share a state pab of the form ([T]) [See 
Fig. [TJa) for reference]. Bob (the sender) wishes to teleport 
(20) an unknown state \i//q) = n\e) + y\g) of an input qubit Q 
to Alice (the receiver). Bob then makes a joint measurement 
on the input qubit Q and the oscillator B, communicates the 
outcome to Alice, who can then implement a correction on her 
qubit A. The final state p^ ut of A can be shown (see below) to 
have a fidelity Ttel = (^qIp^I^q)' averaged over the uniform 
distribution of the input state, given by 

limr xEL =f(l + |r|). (5) 

Notice that for any |r| > the average fidelity Tj EL exceeds 
the classical benchmark achievable by measure-and-prepare 
schemes, T^ L = 2/3 [47]. We can then define a payoff 

*V = (1 - rt)~ l max{0, <F XEL - T^} (6) 

quantifying the better-than-classical performance of the tele- 
portation protocol. It is immediate to see that in the limit of 
large displacements 

P TEL = N(p A B) = 2\r\. (7) 

This shows that entanglement, in the form of negativity, is 
clearly the resource for this protocol [see Fig.[3ja)]. 

We now provide the explicit steps of the protocol and prove 
the result announced in Eq. ([5]), which can be formalized as 
follows: For any e > 0, the state (p} can be used for tele- 
portation of a generic qubit state \i// 1 q) from Bob to Alice, with 

average fidelity Tj Eh > (2/3)(l + |r|)-6 and success probability 

P > 1 - 6. 

We first choose a cutoff integer N and define e N = Yjn+i 
where the initial state of the oscillator B in Eq. ([T]) has been 
written generically as p° B = Ym SnWuX^nV We can find a lower 
bound to the fidelity by just assuming that, with probability 
€n, the protocol fails. That is, we may use p' AB = (1 - ^v)p^ + 
€nPa B where p^ is the truncated version of pab and p^ B a 
state that yields zero teleportation fidelity. From now on we 
assume that Alice and Bob share the truncated resource p^ B \ 
and the derived results will hold with probability P > 1 - e N . 
One can always choose N large enough so that < e. 

The initial state of the complete ABQ system is Pa B q = 
I^qX^qI ®P^5 • The state p^ B can be expanded in terms of the 
qubit basis and the 2(N + 1) states {D(j3)\if/ n ), D 1 (fiM n )} n =o,...,N- 
If we choose large enough, these states will be effec- 
tively orthonormal. Then, Bob performs a hybrid Bell-state 
measurement corresponding to the following 4(N + 1) or- 
thonormal vectors, |0±> = ^(\e)D(fi)\^m) ± \g)D^(J3Mm)l 
and |£> = j=(\e)D^Mm) ± \g)D{J5)Wm)\ These mea- 
surements can be done as follows: Bob performs the dis- 
entangling operation - D\cr^P) on the input qubit Q 
and his oscillator B, then measures the qubit Q in the basis 
{\e) ± and the oscillator in the effectively orthogonal ba- 
sis {\i/f m ) 9 D(2P) ± D\2/3)Mm)}^ =0 US- Then, upon receiving 



4 




FIG. 3: (Color online). Payoffs for hybrid teleportation (a), and for 
hybrid remote state preparation with unitary corrections (b) and op- 
timized with non-unitary corrections (c), plotted versus the correla- 
tions in the shared states, measured by negativity in panels (a),(c) and 
by square root of geometric discord in panel (b); see text for details 
on the boundary curves in (b). The plotted quantities are dimension- 
less and correspond to \J3\ — > oo. 



two classical bits from Bob (as all values of m give the same 
results), Alice first performs a local phase rotation \e) — > j\e) 
on qubit A; further, she may or may not perform the correc- 
tions \e) <-> \g) and \e) — > -\e), depending on Bob's outcomes. 
After tedious calculations, we find that there are only two pos- 
sible (unnormalized) states that Alice gets, 



Clearly the ideal resource for this protocol is obtained for 
p = \r\ = 1/2 and p° B = |0X0|, where |0) is the ground state 
of the oscillator. A typical realistic deviation is given by the 
oscillator being initially in a mixed state. For simplicity, let 
us focus on p° B being in general a Gaussian thermal state l50l 
P°b ~ 2« s n \n){n\, where s n = n n /(l + n) n+l , with n the mean 
number of thermal excitations. This choice allows us to com- 
pute the fidelity analytically (although different choices for p° B 
give rise to qualitatively similar results as verifiable numeri- 
cally), yielding 



lim Trs 



(11) 



In this case, the fidelity is limited by the initial purity of the 
oscillator, akin to the geometric discord in Eq. ([3]). For remote 
state preparation, the classical threshold corresponds to Bob 
preparing a completely random guess (unnormalized) state 
ED El P B l = QPW\ + l"M)/2, yielding 

ni = f d^<v>r > „r \ • (12) 



pT = PM 2 \eXe\ + (1 - p)\y\ 2 \gXg\ + \rW\e\g\ + \r\rfy\gM, Defining again the payoff as 
p°f = (1 - p)\rj\ 2 \eXe\ + p\y\ 2 \gXg\ + \r\rjf\eXg\ + Wy\g\eY 
The input-output fidelity is then given by 

r XEL = J] ^q\pT^q) = M 4 + Irl 4 + 4|r|M 2 | r l 2 . (8) 



(l-O max{0,r RSP -r R c s 



(13) 



, we then find in the limit of large displacements 



L[D G ( PAB )] < P RSP = m3x{0 9 



<yp G (PAB). 



Averaging n, y over the Bloch spere, we obtain Eq. ([5]). 



B. Hybrid remote state preparation 

We now describe a different protocol, see Fig. [TJb) for ref- 
erence. Alice and Bob share again a state pab of the form 
([T]). Without loss of generality, we can assume r = \r\ (this is 
true up to a local unitary on A). Alice (the preparer) wishes to 
remotely prepare l2T1l Bob's oscillator in the target superposi- 
tion state, known to Alice, 



'f) = (10) + e^l-^/^l+cos^e-W), (9) 



for some phase (p. For this purpose, Alice measures qubit A 
in the basis |+a> = (\e) ± e l(p \g))/ V2 and classically communi- 
cates the one-bit outcome "±", obtained with probability P+, 
to Bob. If Alice obtains "+", Bob does nothing, otherwise he 
applies a 7r-phase shift ®# to B, which ideally (for large \(3\) 
maps \/3) — > \P) and -\/3) — > and can be implemented 

in phase space via a combination of displacements and photon 
subtraction |49l . The qumode B after the correction can be in 
one of two possible unnormalized states, p° ut+ = {+a\Pab\+a) 
and p° ut_ = ®b(-a\Pab\-a)®1- The fidelity of the protocol, 
averaged over the distribution of the target state, is then 



(14) 

For a given Dg(Pab), the payoff in the performance of remote 
state preparation can never exceed V^g [which is reached for 
/j° B = 1, solid line in Fig.pTb)], and admits a lower bound as 
well, L[D G ] = max{0, (3lf G - 1)/(1 + D G )} [dashed curve in 
Fig.[5Jb)]. The latter is tight for p = \r\ = 1 /2, when the shared 
states are maximally entangled, yet the remote state prepa- 
ration succeeds with null or limited payoff. A nonzero pay- 
off implies necessarily a nonzero D G . Thus, despite its clear 
shortcomings l25l . the geometric discord might still seem to 
capture the operative performance of qubit-to-oscillator re- 
mote state preparation, in analogy with the case of two-qubit 
resources l24l . 

A remark is in order. In the previously described hybrid 



)K al ). 



(10) 



teleportation (Sec. |III A| ), Bob can optimize his measurement 
strategy to compensate for the initial mixedness of the oscil- 
lator, which therefore does not affect the achievable fidelity. 
In contrast, in remote state preparation the measurement of 
Alice's qubit necessarily leaves Bob's oscillator in a mixed 
state, if its initial state p° B was mixed. As the protocol aims 
at preparing a pure state, it is no surprise that the resource 
mixedness enters in, and degrades, the figure of merit. We can 
address this limitation by extending the conventional remote 
state preparation primitive I2H . to allow Bob to perform gen- 
eral completely positive maps rather than just unitary correc- 
tions, aiming to improve the fidelity with the goal state [51]. 
In our hybrid case, we observe that if Bob performs the opera- 
tion which 'digitalizes' the resource state as discussed above, 
then the mixedness of the oscillator is effectively bypassed 
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and we obtain a payoff = 2\r\, equal to the hybrid tele- 
portation payoff and coinciding with the negativity of pab [see 
Fig.[3jc)]. We now present details of the procedure. 

We start with the usual resource state pab given by Eq. ([I]). 
Alice wants to prepare remotely a state of the form \ifs g B ° al ) = 
-^(|£> + e~ l(p \g)), where \e) and \g) are two orthogonal states 
in Bob's Hilbert space. Alice first rotates her basis such that 
r —> \r\, then performs a projective measurement on her part 
of the system, using the basis \±<p) - ± e~ ltp \g))- The 

reduced density matrix of Bob after the measurement (up to a 
normalization constant) is given by 

(±<p\p AB \±<p) = P D(fi)p° B D^) + (1 - p)rt(0)p° B D(fi)+ 

± \r\e*D(P)p° B D(P) ± \r\e-*rt(P)p° B rt(P). 

(15) 

In the limit of large displacements, Bob can perform the 
following local (non-unitary) operation, with probability ap- 
proaching unity (Note: it is the same 'digitalizing' operation 
detailed in Appendix [C]): 

D(fi)p° B D*(0) -> \eXe\, 
rf(P)p° B D(P) -> 11X11, (16) 
D(fi)p° B D(J3) -> \e\g\. 

The state of Bob's mode then becomes 

p' B = p\m + (1 - P)\m\ ± \A^\e\g\ ± \A&\m> (17) 

which is now properly normalized. Finally, Bob can remove 
the ± signs with a further local unitary correction, obtaining 
eventually the state 

pb = v\m + (i - v)\m\ + \r\**\m\ + vmmv as) 

The fidelity between this output state and the target state is 

n£ = {*7 l \PBW7 l ) = \ + \ri (19) 

which does not depend on the angle and is hence equal to the 
average fidelity T^> • Recalling that for remote state prepa- 
ration T£s P —> \, we finally recover that the payoff, for the 
optimized scheme incorporating a non-unitary local correc- 
tion on Bob's side, becomes - 2\r\ in the limit of large 
displacements, as anticipated above. 

This suggests that the link between remote state prepara- 
tion and measures of discord, highlighted for two-qubit sys- 
tems (23] [24), might be due to a non-optimized version of the 
protocol used. In the case of the qubit-oscillator resources de- 
scribed here, such limitation may be relevant where the 'digi- 
talizing' operation of Eq. ( fT6] ) is experimentally challenging to 
realize, so that one is constrained to unitary corrections only. 
However, the resulting connection between geometric discord 
and protocol performance would be solely due to technolog- 
ical limitations, and not to fundamental quantum-mechanical 
principles. In fact, we showed how properly accounting for 
an extra freedom to correct for the resource mixedness allows 



the protocol to reach a performance only dependent on the 
amount of shared entanglement. 

It will be interesting to test how this result is modified once 
additional decoherence sources in the implementation of the 
hybrid protocols are considered, to see whether the link with 
entanglement will persist also for remote state preparation, or 
some form of discord would emerge as essential operational 
ingredient. 



IV. CONCLUDING REMARKS 

We conclude that, broadly speaking, various indicators of 
"quantumness" in composite systems El [8) can be justified 
by operational interpretations related to their role in differ- 
ent tasks. To wit, entanglement is a resource for teleporta- 
tion (20||47l and superdense coding (52) (among all), geomet- 
ric discord is such for a non-globally-optimized remote state 
preparation scheme l24li (as shown here), entropic discord is 
interpreted through quantum state merging (53) , and relative 
entropy of discord (45] |46) is the cost of entanglement distri- 
bution via separable states (54). The list is likely to grow in 
the future, although care is needed to ensure it does not derails 
away from physical grounds. 

On the practical side, we demonstrated that robust quantum 
correlations in the form of entropic discord and entanglement 
can be engineered in coupled qubit-oscillator systems no mat- 
ter the temperature of the oscillator, provided the interaction 
generates enough displacement. This indicates that, e.g. in 
systems of nano/optomechanical oscillators (55) coupled to a 
two-level probe, there is no need to cool the oscillator down 
to its ground state, in order for quantum communication to be 
achieved efficiently. Although this may require generalized 
(non-unitary) operations that can be challenging to realize 
with current technology, the question of implementing those 
operations in a realistic experimental setup can be a source of 
further interesting research outside the scopes of this paper. 
This may in the future relax the need for ground state cooling 
in favor of generalized operations. In fact, while ground state 
cooling has attracted experimental efforts for a long time, and 
is hence well developed, our findings support the view that the 
experimental realization of generalized quantum operations, a 
comparatively young subject, also may deserve consideration 
and could have a substantial practical impact for the realistic 
implementation of quantum technologies. 
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Appendix A: Negativity of the qubit-oscillator states p A B 

Here we derive Eq. ([2]), which can be formalized as follows. 

Theorem 1. For any e > 0, it is possible to find a fi e (D such 
that 

N(p AB )>2\r\-e, (Al) 

where N is the negativity /fTTT/ (twice the modulus of the sum 
of the negative eigenvalues ofp A A B ). 

Proof The partial transpose of the state ([l}, with respect to 
the qubit, is 

pYb = P\ e M ® D(J3)p° B rt(J3) + (1 - p)\gXg\ rt(P)p%D(P)+ 
+ r\gXe\ D(J3)p° B D(/3) + r>X#l ® D^(fi)p° B D^(fi). (A2) 
Let us consider the test states 

|0 m ) = ^= (k)^^)!^) - e^>D(0)^ m >) , (A3) 

where are the eigenvectors of the initial oscillator state 
P°b ~ ^ml^mX^ml, and is the complex argument of r, that 
is, e^ = j^j . The expectation value of the partial transpose on 
the test states is 

((f>m\p T A A B \(f>m) = ~\r\s m + ^ m \D(2/3)p B D\2/3Mm) 

+ ^<^ m |Z) t (2^Z)(2^)| < A m ). (A4) 

Now, let us make use of the following Lemma (proven below): 
Lemma 2. Given any oscillator density matrix p and any os- 
cillator pure state if/, one has 

lim (i/f\D(a)pD\a)\i/f) = 0. (A5) 

\a\^>oo 

Let us then pick a cutoff N such that Yum Pm > ^ ~ e', where 
e' > 0, and, by using the property ( |A5| ), let us choose J3 such 
that for any m < N one has (ifs m \D(2/3)p B D\2/3)\il/ m ) < e f and 
also (i// m \D\2/3)p B D(2/3)\i// m ) < e' . Then, from Eq. ([A4) it 
follows that 

N 

Yj^AB^rn) < ~\A + <? <\A + N/2). 

m 

By choosing e' = (e/2)(\r\ + N/2)~\ and using the fact that 
the test states are orthonormal, we obtain Eq. ( |A1| ). ■ 

Proof of Lemma 2. Given any density operator p, one has 

lim x P (&) = 0, 

\a\^>oo 

Xp(a) = tr [pD(a)] being its characteristic function. This is a 
consequence of the fact that 

tr [p 2 ]=n~ l f d 2 * p (rf<l. 



In particular, by taking a projector p = \(fX(f\, we have 
lim (<p\D(a)\(p) = 0. 

\a\^>oo 

Then, by substituting \ip) with ± \ip 2 ), |<£>i> ± /|<£>2) in the 
above result, it is easy to show that, for any pair of vectors 
Vpi )> 1^2), one has 

Iim<^i|D(ar)|^ 2 > = 0. (A6) 

\a\^oo 

To prove ( |A5| ), we decompose p = # n |^nX^nl> and write 

(^IZ)^^)^) = J] *J MD(<*)|^>| 2 . (A7) 

We now choose N such that Yn+i Qn ^ 6 /2, and thanks to 
Eq. \A6) we can choose a such that \(if/\D(a)\(p n )\ 2 < e/2 for 
any n < N. It follows that 

2^i<^(«)i^>i 2 

N 00 

= Y j q n mD(a)\<p n )\ 2 + £ q n \(i{,\D(a)\<p n )\ 2 

n n=N+l 

N 00 00 

^\Yu qn+ Yu qn -\Yu qn + \ = e (A8) 

n N+l n 

Then, given any e > 0, we have found a such that 

(i/f\D(a)pD\a)\i/r) < 6, (A9) 
which proves Eq. ( |A5| ). ■ 

Appendix B: Geometric discord of the qubit-oscillator states 

Pab 

1. Geometric Discord for 2 ® 00 systems 

Given a quantum state p^g of a (D 2 (D J system, with 
A being a qubit and 5 being an arbitrary (finite or infi- 
nite) d-dimensional system, the normalized geometric discord 
£>g(Pab) is defined as d El 

£>g(Pab) = 2M\\p AB - U A (p AB )\\ 2 2 , 

ru 

where the infimum is over all von Neumann measurements 
n A = {II*} on the qubit A, with U A (p AB ) = 2^(11* 

WpAB^ k A t B ), and ||M|| 2 = V Tr ( MMt ) is the Hilbert- 
Schmidt norm. 

Noting that any projective measurement induces a dephas- 
ing of the qubit on some orthonormal basis {\c), i.e., 
n A (p AB ) = (p AB + u APAB ul)/2, with U A = \c\c\ - Ic^X^I 
a "root-of-unity" operation on A [56-58], we can recast 
the problem into an optimization over such local unitaries, 
®g(Pab) = 2M Ua \\[p AB - (U A l B )p AB (U\ t B )-\l2\\ 2 r In- 
troducing a generic projector P A = \c\c\ = (e • cr)/2 on the 
qubit, where a = (1, (?) is a four- vector of Pauli matrices, and 
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e = (1, e) with e e R 3 a unit vector, the geometric discord can 
be expressed as follows (57), 

£>g(Pab) = inf (4tr[p^(P>l fi )-p Afi (P>l fi )p Afi (P>l fi )]} . 

(Bl) 

This minimization can be solved in closed form. We define 
a 'partial Fano representation' |59l of p AB by expanding only 
the qubit A in the Bloch basis, 

Pab = ~(v-cr), 

with the four- vector of operators v = tr A (p AB (r) = (vo,v), 
where vo = Pb is the reduced density matrix of B. Like in 
relativity theory we use Greek indices to indicate the com- 
ponents 0, 1, 2, 3 and Roman indices to indicate the 'spatial' 
components 1,2,3 only. 

Let us evaluate each term appearing inside the minimization 



in Eq. (Bl ). For the first, we have 

4tr[pl^] = itr[(vcr) 2 (e-cr)] 

= ^tr fi [v^Vy] e n M A [cT-^o-yCr^] , 

where sum over the repeated indices is understood as in Ein- 
stein's convention. To evaluate the second term, we begin by 
noting that, if P^ = (e • a) = \c\c\, then 

M\PabPaPabPa\ = tv B [tTA \pab\c\c\p ab \c){c\]] 
= M B [(c\p AB \c)(c\p AB \c)] = 
= tr B [tT A [pP AB ]tT A [pP AB ]]. 

Then we have 

4tr [p AB P\p AB P A ] 
= ^tr B [tr A [(v • <r)(e • &)] tr A [(v • <r)(e • <r)]] = 

= itr fi [v^Vy] e^ttA [cr^] tr A [cr y o- T ] = 
= tr s [v^Vy] e^e v . 

For the last equality we have used the fact that tr A [cr^cr^j = 
26 m . We now define the 4 x 4 matrix 

Spn, = M B [v^Vy] . 
Note that the matrix S is symmetric, due to the cyclic in- 



variance of the trace. Then Eq. (Bl) can be rewritten as 
£>g(Pab) = inf e T, with 

T = -S llv e J1 M A [(T^crycr^] - S^e^ey = 

= ^Sije^u [cTiCTjcr^ - Sijetej, (B2) 

where we have used the fact that the contributions with p = 
or y = vanish (this is easy to check in the above formula 



noting that eo = 1). To simplify further the above expression, 
we evaluate explicitly the term 

erjM A [cTiCTjcr^ = 25^ + 2ie ijk e k . 

Since S is symmetric, we have 5^-6^ = 0, and ( |B2| ) simplifies 
to 

T = Su -Sijeiej. 

Defining S = (S)ij as the 3 x 3 'spatial' sub-block of <S, we 
can write 

T = M[S]-fSe. 

Since \e\ = 1, the minimization of such expression yields 
Eq. dB5l, 



£>g(Pab) = HS ) - A max (S ) , with S = M B [W T ] . (B3) 

Recall that S = M B [vv T ], with v = tr A [p AB &]). 

Equation ( |B3[ ) encompasses the known formulas for two- 
qubit C21 I37B and qubit-qudit states El W\ El, but is 
valid as well for states of a qubit and a qumode, for which 
d = oo. In the latter case, it can be convenient to adopt 
a hybrid Hilbert-space/phase-space picture to describe qubit- 
oscillator states p AB (TS1 El). Let b (b^) denote the annihi- 
lation (creation) operator for the qumode, with [b, M] = 1, 
and let D(J3) = exp(j3b^ - J3*b) denote the corresponding 
Weyl displacement operator, with J3 e (D. We can define 
the characteristic vector associated to the state p AB as x(P) = 
(Xo(P)M) = lr\pAB(TD(P)} = M B [vD(J])l where the zero-th 
component Xo(fi) = ^b[PbD(J3)] is the conventional charac- 
teristic function of the oscillator [50], describing its marginal 
state in phase space. The matrix S appearing in Eq. (B3) 
is then S = n~ l J &f3x(P))t(F>)- Similar representations can 
be provided by employing e.g. the Wigner distribution to de- 
scribe the oscillator flIfll4T1|gTI 



2. Explicit calculation for the states p A B 

To calculate Dg(Pab), we need the spatial components of 
the vector v (see previous subsection). These are 

V! = rD(P)p B D(p) + r*rt(P)p B rt(fi), (B4) 
v 2 = irD(J3)p° B D(J3) - ir*rt(JS)p° B rt(JS), (B5) 
v 3 = P D(J3)p° B rt(J3) - (1 - p)rt(fi)p° B D(fi). (B6) 

The matrix S can now be calculated according to = 
us [v/Vy-J. We are interested in the regime of large displace- 
ments. Suppose that is large compared to the phase space 
extension of the initial oscillator state p° B . Then the matrix S 
converges to 



( 2\r\ 2 







2\r\ 2 







2/7(1-/7) + ! 



tr,[^) 2 ]. 



(B7) 
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From the above, as well as the condition \r\ 2 < p(l -p), we can 
directly see that A max (S) = [2p(l - p) + l]tr fi [(p° B ) 2 ]. Hence, 
the geometric discord of the state pab, for the case of large dis- 
placements, is given byEq. 3 \,D G (p AB ) 4|r| 2 tr B (p° B ) 2 . 

\6\— »oo L J 

Note that Dg(Pab) < tr# [(P^) 2 ]* i- e -' me geometric discord is 
smaller than the purity of the initial state of the oscillator. 

To see that Eq. ( |B7| ) is correct in the limit of large dis- 
placements, we can proceed as follows. When we calcu- 
late tr B [v/v 7 J, we get terms of the form tr B \pb£>\F>2PbD?>Da\, 
where each one of the D/s can be either D(J3) or D^(J3). In 
the limit of large \/3\, however, only those that evaluate to 
tr B [(p#) 2 ] survive. All the others can be shown to be negligi- 
ble by using the fact that \tr B [p B D(a)]\ -> for \a\ » \/3\ (see 
also Lemma 2 in the previous Appendix). For example, us- 
ing Glauber's ^-representation we have Pb = f ^aP{a)\a\a\, 
and we can see that (denoting D(J3) = £>), for instance, 

\tr B [p B (rt) 2 p B D 2 ]\< J d 2 ad 2 a f \P{a)P(a')\\{a'\D 2 \a)\ 2 = 

= J d 2 ad 2 a' \P{a)P{a')\z-W +a - a '\ 2 



0, 



(B8) 



and so on. The above and all the other terms can be evaluated 
explicitly e.g. in case p° B is a Gaussian state, but this restriction 
is not crucial for the validity of Eq. ([3]). 



Appendix C: Digitalization of the qubit-oscillator states p A B 

Theorem 3. For any e > 0, it is possible to find a/3 e (D and a 
local operation on Bob, Yjj OjPabO^, with 0^0 7 = 1, such 
that the state pab [Eq. is converted into the 'digitalized' 
two-qubit state 



PAl = p\eXe\®\eXe\ + (l-p)\gXg\® 
+ r\eXg\®\eXg\ + r*\gXe\®\gXe\, 



with probability P > 1 - e and fidelity T > 1 - e. 

Proof. Recall p° B = s n \ifr n X^ n \- Let us fix e N = 
Yi7=N+i s n'i clearly lim^co £n = 0. Let us now con- 
sider the two subspaces *Hi2 spanned by the bases T>\ = 
{D08)|^ri),...,D08)|^)}andS2 = {D(tf\i/r 1 ),... 9 D<fi?\i/r N )}. 
The bases £>i and &2 may overlap with each other; neverthe- 
less, in the limit —> 00, = {Si,^} becomes an orthonor- 
mal set. Let us then choose \/3\ large enough such that this 
orthonormality is verified for all practical purposes. We de- 
note by *Wi2 = *H\ 0*7^2, the subspace spanned by the basis < B. 
Let us consider the local operation on Bob, corresponding to 

Oj = \eXi/sj\rt(P) + VgWiWX 

where j = 0, . . . , N. To see that those operators are quasi- 
complete, we evaluate 

N N 

2 0)0 j = D(fi) WjMj&QS) + D*(fi) \Wj\D(fi) 



P12, 



(C2) 



where *Pn is the projection on the subspace *Hu. Let us in 
general denote the complete operation with {0^ =v ± }. It is 
irrelevant how specifically we complete the operation on the 
remainder of Bob's Hilbert space; we can just assume that the 
fidelity of the output state is zero if we go outside *Hi2 (this 
happens with probability < Now we can easily see that 



7=0 7=0 



lim V O 



so that the output of the complete operation is some state 

Pab = (1 - € n)PaI + 

where p^ B is some state orthogonal to p^J. It is clear that, if 
we choose small enough, this state p^ has fidelity T > 



( c 1 ) l- 6 with the target state p^J of Eq. ( C 1 ) . 
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